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o 

We present numerical calculations of quantum transport in perfect octagonal approximants. 
These calculations include a Boltzmann (intra-band) contribution and a non-Boltzmann 
(inter-band) contribution. When the unit cell size of the approximant increases, the magnitude 
of Boltzmann terms decreases, whereas the magnitude of non-Boltzmann terms increases. It 
shows that, in large approximants, the non-Boltzmann contributions should dominate the 
transport properties of electrons. This confirms the break-down of the Bloch-Boltzmann the- 
ory to understand the transport properties in approximants with very large unit cells, and 
then in quasicrystals, as found in actual Al-based approximants. 

£ 

O 

O . 1. Introduction 



Since the discovery of Shechtman et al. [1| many experimental investigations have 
indicated that the conduction properties of several stable quasicrystals (AlCuFe, 
AlPdMn, AlPdRe...) are at the opposite of those of good metals @, S|. It appears 
also that the medium range order, over one or a few nanometres, is the real length 
O'N scale that determines conductivity. In particular, the role of transition elements 

enhancing localisation has been often studied [I -11]. There is now strong evidence 



that these non standard properties result from a new t ype of break-down of the 



semi-classical Bloch-Boltzmann theory of conduction (l2j, llj] . On the other hand 



the specific role of long range quasiperiodic order on transport pro 



perties is still 



an open question in spite of a large number of studies (Refs. [1, 14-23j] and Refs. 
therein). In this paper, we study "how electrons propagate" in approximants of 
the octagonal tiling. This tiling is one of the well-known quasiperiodic tilings that 
have been used to understand the influence of quasiperiodicity on electronic trans- 
port [Il4l3,[ll421]. The main objective is to show that non standard conduction 



properties result from purely quantum effects due to quasiperiodicity that cannot 
be interpreted through the semi-classical theory of transport. 
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Figure 1. Vertex configurations of the octagonal tiling and of its approximants k > 1. 



2. The octagonal tiling and its approximants 



The octagonal, or Ammann-Beenker, tiling |24l . |25I | is a quasiperiodic tiling analo- 
gous to the notorious Penrose tiling. It has eightfold Cs v point symmetry, inflation 
symmetry with ratio A = 1 + y/2, quasiperiodic translation symmetry, and Bragg 
peaks indexed by Z 4 . Those symmetries were observed in quasicrystalline phases 
of CrNiSi and AlMnSi alloys 0]. 

A sequence of periodic approximants Xq,Xx, . . . , X k , . . . can be generate d by 
the inflation mapping M, operating as an invertible integer matrix on Z 4 27]. 



In the physical Euclidean plane E, and in the simplest version where the atoms 
occupy the tiling vertexes, all the (atomic) sites have integer coordinates in the 
basis (ei, e2, e3, 64), where (ei,e2) forms an orthogonal basis of E and (63,64) is 
another orthogonal basis related to the former by 45° rotation. All the tiles are 
congruent to either a square or a rhombus with 45° acute angle. 

For k = 0, 1,2, . . ., the /c th approximant is periodic, invariant under the trans- 
lation lattice C k generated by \ k e\ and A fc e2 (the square lattice [ei,e2] dilated 
by a factor L k = \ k ). In terms of the Pell, or "octonacci", number sequence 
P = {Pk} = {0, 1, 2, 5, 12, . . .}, defined by P = 0, P x = 1, P k+1 = 2P k + 
the unit cell contains N k = i^Jfc+l + P2k vertexes. Those vertexes are given by the 
following construction: \ k (W k H [63,64]), i.e. select and expand sites of the square 
lattice [ea, 64], eventually bring them back into the square fundamental domain by 
C k translations. The selection window W k is a semi regular octagon obtained as a 
projection of the 4D unit cube. Up to a sign (-l) fc , W fc is the span (all combinations 
with coordinates between and 1) of the four vectors 



Pk{e^ ~ ez), -P k (e 3 + e 4 ), (P k + Pfc_i)e 3 , (P k + Pfc_i)e 4 . 



Ref. [27j provides a slightly more efficient numerical method, avoiding useless 
scans, by partitioning the octagon into six parallelogram sub-domains; then each 
selection reduces to specifying ranges for the integer coordinates of a suitable lat- 
tice, followed by a global affine transformation. 

The first approximant Xq is a simple square tiling. In all subsequent approxi- 
mants k > 1, the configurations around vertexes are the same as in the octagonal 
quasiperiodic tiling, displayed in Fig. [TJ Only the approximants X\ contain incom- 
plete subsets of these six local patterns, depending on the phase. 

In units of edge length a = ||ej||, the shortest distances between vertexes are 
l/v2j 1, V%, corresponding to the short rhombus diagonal, the tile edge and the 
square diagonal respectively. 
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Figure 2. (colour online) Total density of states in perfect octagonal approximants k = 3 and k = 5. DOS 
is symmetric w.r.t. E = [20ll . 
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Figure 3. (colour online) Average Boltzmann velocity Vb in perfect octagonal approximants k = 3, k = 4 
and A; = 5. 



3. Electronic structure 



We study a pure hopping Hamiltonian 



H 



(i) 



where i indexes s orbitals located on all vertexes and 7 = 1 eV is the strength of the 
hopping between orbitals. are the nearest-neighbours at tile edge distance a 
(figure H| . The properties of this model depend only on topology of the tiling. The 
total density of states (total DOS), n(E), is computed by a recursion method 22] 
in a large super-cell containing more than about 10 million orbitals. DOS is shown 
in figure [21 As already shown by E. S. Zijlstra [2(|, DOS is spiky at the centre of 
the band (\E\ < 2) and smooth near the band edges (\E\ > 2). 

The Boltzmann velocity (intra-band velocity) Vb oc dE Qq^ (q- wave vector) is 
shown figure El When the size of the approximant increases, Vb decreases as ex- 
pected from band scaling analysis In order to quantify this localisation phe- 
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nomenon, we compute the participation ratio of each eigenstate \tp) defined by 

N _j 

= (Nj2\m\ A ) , (2) 

where i indexes orbitals in a unit cell and N{= N^) is the number of atoms in 
this unit cell. For completely delocalised eigenstates p is equal to 1. On the other 
hand, states localised on one site have a small p value: p = 1/N. The average 
participation ratio p{E) of all states with energy close to E are presented on figure 
2J This figure shows clearly a stronger localisation of electronic states for larger 
approximants (larger k values). 



4. Quantum diffusion 

In the framework of Kubo-Greenwood approach for calculation of the conductivity, 
the central quantities are the velocity correlation function of states of energy E at 
time t, 

C(E,t) = (v x (t)V x ($) + V x (W x {t)) e = 2Re(v x (t)V x (0)) E , (3) 

and the average square spreading of states of energy E at time t along the x 
direction, 



AX 2 (E,t) = ((X(i)-X(O) 



(4) 



In ([3]) , Re A is the real part of A and V x (t) is the Heisenberg representation of the 
velocity operator V x along x direction at time t: 



in 



X,H 



(5) 



C(E,t) is related to quantum diffusion by: 

±(AX 2 (E,t)) = J*C(E,t')dt'. 



(6) 
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Figure 5. (colour online) Velocity correlation function C(E,t) versus time t in perfect octagonal approx- 
imants k = 3, k = 4 and k = 5. For E = 0.5 eV. 



At zero temperature, the static conductivity is given by the Einstein formula, 

a = e 2 n(E F )D(E F ) , (7) 

where E F is the Fermi energy and D, the diffusivity. In a perfect tiling at temper- 
ature T = K, 

D [E) = AX2 ( E '*) ( as t -> +00). (8) 

Once the band structure is calculated from the tight-binding Hamiltonian ([TJ the 
velocity correlation function can be computed exactly in the basis of Bloch states 



12, l28j. Relation ([6]) shows that an anomalous behaviour of C(E, t) also implies an 
anomalous behaviour of the quantum diffusion. In crystalline approximants, at the 
long time, when the spreading of states, AX = V AX 2 , reaches the size of a unit 
cell, the propagation should become ballistic. This means that the time averages 
at long times are C(E,t) ~ 2Vg and AX 2 (E,t) ~ Vgt 2 , with V B the Boltzmann 
velocity along the x direction. But for smaller time values, when AX < the 
spreading of states depends only on the tiling structure within a unit cell, which is 
closely related to quasiperiodicity. Therefore, an analysis of electron propagation 
in approximants on time range before the ballistic regime gives a good indication 
on how electrons propagate in quasiperiodic structure. 

As shown on figure O the velocity correlation function in approximants can take 
negative values at some times, indicating the presence of back-scattering effects in 
electron propagation. The decay of the Boltzmann term when k increases suggests 
that, for larger approximants, back-scattering (C(E,t) < 0) should be obtained 
even at some large time values. In that case, conductivity of the tiling will be 
strongly affected by this back-scattering phenomenon (equation [8]) . 

The average square spreading is the sum of two terms (figure [6|) : 

AX 2 (E,t) = Vit 2 + AX 2 B (E,t). (9) 

The first term is the ballistic (intra-band) contribution at the energy E. The semi- 
classical model of Bloch-Boltzmann theory is equivalent to taking into account 
only this first term. The second term (inter-band contributions) AX^ B (E,t) is a 
non-ballistic (non-Boltzmann) contribution. It is due to the non-diagonal elements 
of the velocity operator and describes the spreading of the wave-function in a cell. 
This term is bounded by L m (E) 2 , which depends on the energy e H3|. From nu- 
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Figure 6. (colour online) Boltzmann (AJfg) and Non-Boltzmann (AX^ B ) square spreading versus time 
t at E = 0.5 eV in perfect octagonal approximants k = 3, k = 4 and k = 5. a is the length of edges of the 
tiling. 
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Figure 7. Non-Boltzmann square spreading AJf^ B versus time t at different energies E in perfect octag- 
onal approximant k = 5. a is the tile edge length. 



merical calculations, it is found that AX^ B (E, t) reaches rapidly its maximum limit 
~ L m (E) 2 (figureEl). Moreover for many energies, numerical calculations show that 
L m {E) is of the order of magnitude of the unit cell size L/% or less (see also Ref. [12] )• 
It suggests that in approximants the medium range order localises partially electron 
states and this phenomenon should lead to anomalous diffusion in quasicrystalline 
phases. For large time (t > L m /Ve) the ballistic contribution dominates but for 
small time (t < Lm/Vs) the non-ballistic contribution could dominate ("Small 
velocity regime" [12]). When the size of the approximants increases, the charac- 
teristic time L m /VB of the crossover between ballistic and non-ballistic behaviour 
increases. 

In the relaxation time approximation 0, EH, 0, H§] , the role of phonons and 
static defects are taken into account by a scattering time r. r decreases when 
temperature T increases and when the number of static defects increases. The 
velocity correlation function C at temperature T and/or with static defects is 
computed from C at T = K in perfect tiling, 



C(E,t) = C(E,t)e 



\t\/r 



(10) 



Here the Fermi-Dirac distribution function is taken equal to its zero temperature 
value. This is valid provided that the electronic properties vary smoothly on the 
thermal energy scale k-gT. The static diffusivity D' is then, dF 



i r+oo 

D'(E f ,t) = - / C(E F ,t)e-^ T dt. 
1 Jo 



(11) 
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Figure 8. (colour online) DifTusivity D' , in relaxation time approximation, versus scattering time r, at 
Ep = 0.5 eV in octagonal approximant k = 5. D 1 = DL + DL B . 



Roughly speaking, the transport properties are almost entirely determined by the 
AX 2 (E,t = t). Therefore, in large approximants for which (r < L m /V&), trans- 
port properties will be governed by non-ballistic (non-metallic) behaviour. Figure 
[5] shows the diffusivity in approximant k = 5, calculated in the relaxation time 
approximation. D' is the sum of a ballistic term, D' B oc t, and a non-ballistic term, 
-D^B- F° r small r value, its behaviour is not ballistic (D 1 ^ D' B ). In larger approx- 
imants (k > 5), one can expect an "insulating like" behaviour (D' decreases when 
r increases) for a realistic range of r values, as that found in real approximants 



5. Conclusion 

To summarise, we have presented quantum diffusion by a pure hopping Hamiltonian 
in approximants of octagonal tiling (k, < 5). When the size of the approximant 
increases the usual Boltzmann term decreases rapidly and new non-Boltzmann 
terms become essential to understand transport properties. These non-Boltzmann 
terms can have "insulating like" behaviour, suggesting that in larger approximants 
(k > 5), "insulating like" states, due to long range quasiperiodic order, could exist. 
Calculations in larger approximants are in progress. 
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